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The kink stability of self-similar solutions of a massless scalar field with circular symmetry in 
2 + 1 gravity is studied, and found that such solutions are unstable against the kink perturbations 
along the sonic line (self-similar horizon). However, when perturbations outside the sonic line are 
considered, and taking the ones along the sonic line as their boundary conditions, we find that non- 
trivial perturbations do not exist. In other words, the consideration of perturbations outside the 
sonic line limits the unstable mode of the perturbations found along the sonic line. As a result, the 
critical solution for the scalar collapse remains critical even after the kink perturbations are taken 
into account. 
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I. INTRODUCTION 

The studies of non-linearity of the Einstein field equations near the threshold of black hole formation reveal very 
rich phenomena [1], which are quite similar to critical phenomena in statistical mechanics and quantum field theory 
t-H , [2]. In particular, by numerically studying the gravitational collapse of a massless scalar field in 3 + 1-dimensional 
< ^ • spherically symmetric spacetimes, Choptuik found that the mass of such formed black holes takes the form [1], 
\Q ■ 

O M BH = C(p)(p-p*y , (1-1) 

in . 

where C(p) is a finite constant with C(p*) ^ 0, and p parameterizes a family of initial data in such a way that when 
""^ p > p* black holes are formed, and when p < p* no black holes are formed. It was shown that, in contrast to C(p), 
O"^ the exponent 7 is universal to all the families of initial data studied, and was numerically determined as 7 ~ 0.37. 

The solution with p — p*, usually called the critical solution, is found also universal. Choptuik's studies were soon 
6JT). generalized to other matter fields [3]. From all the work done so far, the collapse in general falls into two different 
types, depending on whether the black hole mass takes the scaling form (1.1) or not. When it takes the form, the 
corresponding collapse is called Type II collapse, and when it does not it is called Type I collapse. In the type II 
collapse, all the critical solutions found so far have either discrete self-similarity (DSS) or homothetic self-similarity 



(HSS), depending on the matter fields. In the type I collapse, the critical solutions have neither DSS nor HSS. For 
certain matter fields, these two types of collapse can co-exist. A critical solution in both two types has one and only 
one unstable mode. This now is considered as one of the main criteria for a solution to be critical. 

The studies of critical collapse have been mainly numerical so far, and analytical ones are still highly hindered by 
the complexity of the problem, even after imposing some symmetries. Lately, some progress has been achieved in 
the studies of critical collapse of a scalar field in an anti-de Sitter background in 2 + 1-dimcnsional spacetimes both 
numerically [4,5] and analytically [6-9]. This serves as the first analytical model in critical collapse. In particular, 
Garfinkle [6] first found a class of exact solutions to Einstein-scalar field equations, denoted by S[n], and later Garfinkle 
and Gundlach (GG) studied their linear perturbations and found that the solution with n = 2 has only one unstable 
mode [8]. By definition this is a critical solution, and the corresponding exponent 7 in Eq.(l.l) can be read off from 
the expression [10] 
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from which it was found 7 = 4/3, where ki denotes the unstable mode. Although the exponent 7 is close to that 
found numerically by Pretorius and Choptuik [4], 7 ~ 1.2 ± 0.05 (but not to the one of Husain and Olivier, 7 ~ 0.81), 
this solution is different from the numerical critical solution [6]. Using different boundary conditions, Hirschmann, 
Wang and Wu (HWW) found that the solution with n = 4 has only one unstable mode [9] . As first noted by Garfinklc 
[6], this n = 4 solution matches extremely well with the numerical critical solution found by Pretorius and Choptuik 
[4]. However, the corresponding exponent 7 now is given by 7 = = 4, which is significantly different from the 

numerical ones. The boundary conditions used by HWW are [9]: (a) The perturbations must be free of spacctimc 
singularity on the symmetry axis; (b) They are analytical across the self-similarity horizon, as the background solutions 
do; (c) No matter field come out of the already formed trapped region [11]. GG considered only Conditions (a) and 
(b) [8]. 

In this paper we shall study another important issue of critical collapse for a scalar field, the kink stability. The 
kink modes result from the existence of critical characteristic lines (they are also referred to as self-similarity horizons, 
and sonic lines), along which discontinuities of (higher order) derivatives of some physical quantities can be developed 
and propagate. The instability is characterized by the divergence of the discontinuity, and the blow-up may imply 
the formation of shock waves [12]. An example that discontinuities of derivatives can propagate along a sonic line is 
given by the linear perturbation, Sip(r 7 z) — ipi(z)e kT , of the massless scalar field in 2 + 1 gravity, which satisfies the 
following equation [9], 

z(l-z) W + \ [(1 + 2fc) - z(3 + 2k)} Vl ' - ifc^i = f(z), (1.3) 

where a prime denotes the ordinary differentiation with respect to the indicated argument, f(z) is a smooth function 
of z, and z = 1 is the location of the sonic line (cf. Eq.(112) in [9]). From the above equation we can see that it 
is possible for ipi to has discontinuous derivatives only across the line z — 1. In fact, assume that ipi is continuous 
across z = 1 (but not its first-order derivative), we can write it in the form 

Vl {z) = p+(z)H(z - 1) + <pi(z) [1 - H(z - 1)] , (1.4) 

where H (x) denotes the Heavside (step) function, defined as 



Then, we find that 



Vl ' = l p+'H(z-l) + ^'[l-H(z-l)], 

Vl " = <fit"H(z - 1) + <pf [1 - H{z - 1)] + [ipi'Y 5{z - 1), (1.6) 



where S(x) denotes the Dirac delta function, and 



Substituting Eq.(1.6) into Eq.(1.3) and considering the facts 

H m (x) = H(x), [1 - H{x)] m = [1 - H(x)} , 

[1 - H(x)] H(x) = 0, x6(x) = 0, (1.8) 

where m is an integer, one can see that Eq.(1.3) holds also on the horizon z = 1 even when [<pi] ^ 0. This is because 
(1 — z) [</VP S(z — 1) = 0, as long as [</VP is finite. Thus, we have 

[^,J~ = [^i're fer , (1-9) 

where ( ) iZ = d( )/dz. The above expression shows clearly how the discontinuity of the first derivative of the 
perturbation Sip(r, z) propagate along the sonic line z = 1. When Re(k) > the perturbation grows exponentially 
as r — * 00, and is said unstable with respect to the kink perturbation. When Re(k) < the perturbation decays 
exponentially and is said stable. 

Note that if the discontinuity happened on other places, say, z = z ^ 1, clearly Eq.(1.3) would not hold on z = z , 
because now (1 — z) [(p\]~ S(z — z ) ^ 0. This explains why the discontinuities are allowed only along the sonic lines. 
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The above analysis also shows that the kink perturbations are different from the ones considered in [8] and [9] , because 
there it was required that ipi(z) is analytical across z — 1, that is, 



0, (to -1,2,...) (1.10) 



where tpi^ denotes the m-th order derivatives of ip\. Therefore, kink perturbations were excluded in the studies of 
linear perturbations of [8] and [9] . 

Ori and Piran first studied kink stability of self-similar solutions in newtonian gravity [13], and lately Harada 
generalized such a study to the relativistic case and found that the critical self-similar solutions of a perfect fluid 
with the equation of state P — kp are not stable against kink perturbations for k > 0.89, where P and p denote, 
respectively, the pressure and energy density of the fluid [14]. More recently, Harada and Maeda showed that in 
four-dimensional spherically symmetric case the self-similar massless scalar solution found lately by Brady et al [15] 
is also not stable against kink perturbations [16]. 

In this paper, we study the kink stability of the scalar field in 2 + 1 gravity. Instead of assuming that 8<p(t, z) is 
c° across the sonic line, as we did in the above example, following Harada [14], and Harada and Maeda [16], we shall 
assume that S(p(r, z) is c , that is, 5(p(r, z) and its first-order derivative with respect to z are continuous across the 
sonic line, but not its second-order derivative. We shall first show that perturbations obtained along the sonic line 
allow the existence of unstable modes. However, when we consider perturbations outside the sonic line and take the 
ones obtained along the sonic line as their boundary conditions, we find that these conditions together with the ones 
on the symmetry axis do not allow any non-trivial perturbations in the regions outside the sonic line. Therefore, the 
consideration of perturbations in the whole spacetimc limits the unstable mode found along the sonic line. Thus, all 
the self-similar solutions of the massless scalar field are stable against kink perturbations in 2 + 1 gravity. As a result, 
the critical solution for the scalar collapse remains critical, even after the kink perturbations arc taken into account. 

Specifically, the paper is organized as follows: In Sec. II we give a brief review of the self-similar solution, which 
is needed in the studies of linear perturbations in Sec. Ill, in which we first consider the linear perturbations of the 
self-similar solutions along the sonic line, and then the perturbations outside the sonic line. In Sec. IV, we summarize 
the main results obtained in this paper and then present our concluding remarks. 



II. THE EINSTEIN-SCALAR FIELD EQUATIONS 

The general form of metric for a (2 + l)-dimensional spacetime with circular symmetry can be cast in the form, 

als 2 = -2e 2a{u ^ dudv + r 2 {u, v)d8 2 , (2.1) 

where (u, v) is a pair of null coordinates varying in the range (— oo, oo), and 9 is the usual angular coordinate with 
the hypersurfaces 9 — 0, 2tt being identified. = dg is a Killing vector. It should be noted that the form of the 
metric is unchanged under the coordinate transformations, 

u = u(u), v — v(v). (2.2) 

To have circular symmetry, some conditions on the symmetry axis needed to be imposed. In general this is not trivial. 
As a matter of fact, only when the axis is free of spacetime singularity, do we know how to impose these conditions. 
Since in this paper we are mainly interested in gravitational collapse, we shall assume that the axis is regular at the 
beginning of the collapse. In particular, we impose the following conditions: 
(i) There must exist a symmetry axis, which can be expressed as 



X 



0, (2.3) 



as v — > f(u), where we assumed that the axis is located at r(v = f(u),u) = 0. 

(ii) The spacetime near the symmetry axis is locally flat, which can be written as 



^ - 1, (,4> 



as v — > f(u), where ( ) iQ = d( )/dx a . 

The corresponding Einstein-scalar field equations for the metric (2.1) take the form, 
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r tUU - 2<7„r „ = -8irGr(f> 2 ul (2.5) 

r tVV - 2a. v r. v = -8nGr4? v , (2.6) 

r, uv + 2ra tUV = -8nGr(f), u <j) tV , (2.7) 

r.uv = 0, (2.8) 
while the equation of motion for the scalar field is given by 

20,™ + - r (r, + r „£„) = 0. (2.9) 
To study self-similar solutions, we first introduce the dimcnsionless variables, z and t, via the relations 

= -m(^l, (2.10) 



(-«)' \ "0 

where w is a dimensional constant with the dimension of length, and the above relations are assumed to be valid 
only in the region v > 0, u < 0. We will refer to this region as Region / [cf. Fig. 1]. 
Self-similar solutions are given by 

F(t,z)=F ss (z), (2.11) 

where F = {a, s, ip}, and 

r(u,v) = (-u)s(r, z), 

(p(u,v) = cln\-u\ + ip{T,z), (2.12) 

with c being an arbitrary constant. A class of such solutions was first found by Garfinkle [6], which can be written as 
[9] 

= -In | ^—^ |+a , 

r ss (u, u) = (-u) - v, 

<j> ss (u, v) = 2c In Yu 1 ' 2 + e(-u) 1/2 | + (2.13) 

where e = ±1, (Tq and </?J are integration constants, and \ = 8ttGc 2 . As shown in [9], the hypersurface v = for 
the solutions with 1 > % > 1/2 represents a sonic line, and the solutions can be extended across the hypersurface, 
whereby they can be interpreted as representing the gravitational collapse of a scalar field, in which a black hole is 
finally formed. The extension can be realized by introducing two new coordinates u and v via the relations 

u^-i-u) 1 ' 211 , v = v^ 2n , (2.14) 

where n = 1/[2(1 — x)] > 1. In order to have the extension unique, we require that it be analytical across the 
hypersurface v = 0, which, in turn, requires n to be an integer and satisfy the condition, 

"5Ji^J-{£fi. I = <" 5 > 

where I is another integer. For the detail, we refer readers to [9]. In these new coordinates, the metric and the massless 
scalar field are given by 

ds 2 = _2e 2ff -(°< s ) dudv + r ss 2 (u, v)d9 2 , 
a ss (u,v) = a ss (u,v) + ^\n^4n 2 (-to) 2 " -1 j 

= iln{4n 2 |/(^)| 4 *}+a \ 

r ss (u,v) = (~u) 2n -v 2n 

4> ss {u,v) = 2c\n\f(u,v)\+<t>l, (2.16) 

where f(u,v) = v n + e(—u) n . Note that the symmetry axis (the vertical line r = in Fig. 1) is located at v = u, for 
which conditions (2.3) requires (Tq = ^(1 — 4%) ln(2). The corresponding Penrose diagram is given by Fig. 1. 
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FIG. 1. The Penrose diagram for the solutions given by Eq.(2.16). <f>, a is timelike in Region II and null on the hypersurface 
v — 0. In Region I it is spacelike for n = 21 and timelike for n = 21 + 1. The horizontal line r = is singular for n = 21 + 1 but 
not for n = 21. In Region / all the rings of u, v = Constant are trapped, while in Region II none of them is trapped. 



III. LINEAR PERTURBATIONS OF SELF-SIMILAR SOLUTIONS: KINK STABILITY 

In this section, we consider the linear perturbations of the self-similar solutions given by Eq.(2.16). For the sake of 
simplicity, we shall drop all the bars from ct, u and v, so that the background solutions can be written as, 

ds 2 = -2e 2a ^ u ^dudv + r 2 ss {u, v)d6 2 , 

a ss (u,v) = iln (4n 2 |/(u,v)| 4x ) +a], 

r ss (u,v) = (-u) 2n -v 2n , 
0ss(u,v) = 2cln|/(u,w)| + (j>l, 
f(u,v)=v n + e{-u) n . (3.1) 

Let us first divide the spacetime in Fig. 1 into three different regions, fi ± and S, defined, respectively, by fi + = 
{x a : u, v < 0, v > u} 7 fl~ = {x a : u < 0, v > 0, v < \u\} 7 and S = {x a : v = 0}. Then, for any given C 1 function 
f(u, v), we can write it as 

/(«, v) = f+(u, v) [1 - H(v)} + /"(«, v)H(v), (3.2) 
where f^ 1 denote the functions, denned, respectively, in the regions f^. In the present case, we have 

f ± (u,v) = ft s (u,v)+Sf ± (u,v), (3.3) 
where f ss = {<J SS , r ss , <j> ss } denotes the background solutions given by Eq.(3.1), which are analytical across v = 0, 



Since f(u,v) is C 1 , we must have 



lim Sf + (u,v)= lim <5f (u,v) = 5f c (u), 
v— >o- ' t)-»0+ 

lim */+(«,«) = lim $/-(«,«) = «/ c (1) («). (3.5) 

v— >0~ u— »0+ 



Then, wc find 
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/,„(«, v) = /+(«, w) [1 - + /-(u, w)JT(v), 

/, w (u, w) = /+>, v) [1 - + f; vv (u, v)H{v). (3.6) 

Inserting Eqs.(3.3)-(3.6) into Eqs.(2.5)-(2.9) and considering Eq.(1.8), to the first order of 5f, we obtain 

5r iUU - 2 ((Jss^Sr.u + r ss , u 5a, u ) = ~8ttG (2r ss <j) ss , u 5<j) tU + (f) 2 ss , u 5r) , (3.7) 

5r tVV - 2 (a ss , v 5r, v + r SStV 5a^ v ) = -8ttG (2r as <l> aStV 6<j> iV + 4% a , v 6r) , (3.8) 

2 (r ss 5a iUV + (J S sMv5r) = -8irG [r ss (4> SSlU 5^,v + 4>ss,v54>, u ) + <j>ss,u4>ss,v5r] , (3.9) 

5r, uv = 0, (3.10) 

2r ss 5cp iUV + 2<fi SSiUV 5r + r ss ^ u 5<p, v + r SSiV 5cp iU + (fi ss , u 5r !V + (f> ss ,v5r iU = 0, (3-11) 

where the quantities f ss and 5f should be understood as /+ and 5f + in + , and /~ and Sf~ in f2~. 



A. Kink Stability 

Kink stability is the study of the linear perturbations of Eqs.(3.7) - (3.11) along the sonic line v = 0. To solve these 
equations for 5f c (u), following Ori and Piran [13] (See also [14,16]), we impose the following conditions: Assume that 
the perturbations turn on at the moment, say, u = uo, then we require 

(A) the perturbations initially vanish in the interior, 

Sf-{u ,v) =0,ue fT, (3.12) 

(B) the perturbations and their first-order derivatives be continuous everywhere, and in particular across the sonic 
line, 

[<*/]"= 0, [*/,„]" = 0, (v = 0) , (3.13) 

(C) 5^ vv and 5a ± ^ vv be discontinuous across the sonic line, 

5<j>l{u) = [5<j> iVV \- ± 0, 5a" c (u) = [Sa, vv ]- ^ 0, (v = 0) . (3.14) 

From the above we first note that Eq.(3.12) remains true for all u > u . In fact, Sf~(u,v) = are indeed solutions 
of Eqs.(3.7) - (3.11) in O". Then, from Eqs.(3.12) and (3.13) we find 

Sf+(u,0) = 0, */+(«, 0)=0, 

6^(u) = 54>+ vv (u,0), Sa'^u) = 5a+ v (u,0). (3.15) 

Taking the limit v — > 0~ in Eqs.(3.7)-(3.11) and considering the above equation we find that 

[6r, vv }- =6r+ v (u,0) = 0. (3.16) 

On the other hand, taking derivatives of Eqs.(3.11) and (3.9) with respect to v, and then taking the limit v — > CP, 
we obtain 

2r ss (6^) iU + r ss , u 5<f>'t = 0, (3.17) 

(*O lU = -*xG<f> 6a ,„5<%, (3.18) 

along the sonic line v — 0. Substituting Eq.(3.1) into the above equations and then integrating them, we obtain 

A A 



*c(«) = 



(-«) n u y 2 



72 



**<«> - - (3.19) 

where A is an integration constant. Since n > 1, from the above expressions we can see that both 5<p"(u) and 
5a" (u) diverge as tt — > 0~ (or r — > oo), or in other words, the self-similar solutions are not stable against the kink 
perturbations. 
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It should be noted that Sf + (u,v) cannot be zero identically in f2 + , because we already have 5f (u,v) = in Q 
and 

Sf' c '(u)=5f+ v (u,0-)^0. (3.20) 

Then, a natural question rises: Do the perturbations given by Eq.(3.19) match to the ones in region f2 + ? To answer 
this question, in the next subsection we shall consider the linear perturbations of Eqs.(3.7)-(3.11) in region by 
considering Eqs.(3.15), (3.16) and (3.19) as their boundary conditions at v = 0. 



B. Linear Perturbations in Q + 



To study the linear perturbations in it is found convenient to use the dimensionless variables r and z, defined 
by Eq.(2.10). However, they are valid only in region 0~. In region f2 + we define them as 

f = -]*(=*), ~ z = l, (3.21) 

where u, v < in and 

u = -(-u) 2, \ v = -(-v) 2n . (3.22) 

The null coordinates u and v in Eq.(3.22) should be understood as the ones, u and v, defined by Eq.(2.14). In terms 
of f and z, the background solutions (3.1) in £1+ can be written in the form, 



s (z) = 1 - z, 

<ro(~z) = 2 X \n (V^ + z- 1 / 4 ) 



<p {z) = 2cln(l + z 1 / 2 ) +tp\, (3.23) 

with 

a ss (u, v) = cr (5) + 1 In [Wi-uvf 11 - 1 ] , 
r ss (u, v) = (-u) 2n s a (z), 

(j> ss {u, v) = <p Q (z) + cln [(-u) 2n ] . (3.24) 

Again, u and v in Eq.(3.24) should be understood as u and v defined by Eq.(2.14). For detail, we refer readers to 
Eqs.(57)-(59) in [9]. Without causing any confusions, in the following we shall drop the tildes from f and z. Then, 
writing the perturbations as 

Sr = (-u)si(z)e kT , 
6a = <Ji(z)e kT , 

6<t> = ^i{z)e k \ (3.25) 

it can be shown that the linearized perturbations given by Eqs.(3.7)-(3.11) reduce exactly to the ones of (67)-(71) 
of [9], the general solutions of which are Eqs.(110)-(118) for k = 1, and Eqs.(120)-(125) for k ^ 1, given in [9]. In 
particular, s\(z) is given by 

Sl{z) -\^- k + sl Ml, ( ' 

where (3 and s? are the integration constants. 

However, since here we consider the kink stability, the boundary conditions are different from the ones used in 
[8,9]. In particular, in [8,9] it was required that the perturbations be analytical across the surface v = 0, while in the 
present case these conditions should be replaced by Eqs.(3.15), (3.16) and (3.19), which can be written as 
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k ~ 2' 
Sl (z)~0(z 3 ), 

^(z)^-^z 2 + 0(z 3 ), (3.27) 
2u 

as z — > 0. Thus, the instability of the perturbations along the hypersurface v — found in the last subsection is due to 
a single mode, k = 1/2. In addition to the above conditions, we also need to impose some conditions on the symmetry 
axis r — 0, so that the local-flatness conditions (2.3) and (2.4) are satisfied. In terms of /i(z), these conditions are 
exactly the ones given by Eq.(105) in [9], 

*lWL=l 

a i( z )\z=i ~ finite 



(l-,)^M-2^ l( ,) 



- finite, (r = 0). (3.28) 

2=1 



From Eqs.(3.26) and (3.28) we find that sf = —(5, while Eq.(3.27) requires (3 — 0, for which the solutions of <Ji{z) 
and (fii(z) with k = 1/2 are given by [9] 



c 



jl/2(1 + z i/ 2) #i(f) + l 



dz ' 2 </ ' 1 



(3.29) 



fiiz) = c\F Q, 1; 1; zj +c 2 F Q, i; 1; 1 - zj , (3.30) 

where c\ and C2 are two arbitrary constants, and F (a, 6; c; z) denotes the ordinary hypergeometric function with 
F (a, b; c; 0) = 1. From the expression [17], 

F (i § ^iSr 1 {*<- + '» - * (" + 3) - 5 - "} c - =»"■ (3 - 31 » 

we find 

F {\'\' 1 ' Z ) -~Ml-«). (3-32) 
as z — > 1. Then, Eq.(3.28) requires ci = 0. On the other hand, from Eq.(3.31) we also find 

f(^;1;1-*)->~M*), (3-33) 

as z — > 0. Thus, the conditions of Eq.(3.27) yield c 2 = 0. In review of all the above, we find that the boundary 
conditions (3.27) and (3.28) require 

si{z) = ai(z) = <pi{z) = 0. (3.34) 

That is, non-trivial perturbations in + are not allowed by the boundary conditions (3.27) and (3.28). Then, we 
must have 5f"(u) = 0. In other words, the consideration of the perturbations in fl + limits the unstable mode of the 
perturbations along the sonic line v = 0. 

IV. CONCLUSIONS 

In this paper, we have studied the kink stability of the self-similar solutions of a massless scalar field in 2 + 1 gravity, 
and found that perturbations along the sonic line (self-similar horizon) indeed allow the existence of an unstable mode. 
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In the study of kink stability, it is assumed that the spacetime inside the sonic line is not perturbed, that is, 
Sf~(u,v) = identically [13,14,16]. Then, Sf + (u,v) must be non-zero outside the sonic line, in order to have non- 
vanishing perturbations along the sonic line. However, the perturbations outside the sonic line cannot be arbitrary. 
In particular, they have to match to the ones along the sonic line. In addition, they need also satisfy some physi- 
cal/geometrical conditions, such as, the local-flatness conditions on the symmetry axis. A natural question now is: 
After considering all these, does the spectrum of the perturbations obtained along the sonic line still remain the same? 

To answer this question, in Sec. Ill we have studied the perturbations outside the sonic line, by taking the ones 
obtained along the sonic line as their boundary conditions. We have shown explicitly that these conditions, together 
with the ones on the symmetry axis, indeed alter the spectrum of the perturbations along the sonic line, and in 
particular, they limit all the unstable modes. Thus, all the self- similar solutions of the massless scalar field in 2 + 1 
gravity is stable against kink perturbations. As a result, the critical solution for the scalar collapse remains critical 
even after the kink perturbations are taken into account. 

Finally, we note that in the newtonian gravity the spectrum of the perturbations along the sonic line remains the 
same, even after the perturbations outside the sonic line are taken into account [18]. It would be very interesting to 
see if this is still the case in four-dimensional spacetimes in the framework of Einstein's theory of gravity. 
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